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Abstract 

Kalman filters are often used to estimate the state variables of a dynamic system. 
However, in the application of Kalman filters some known signal information is often 
either ignored or dealt with heuristic-ally. For instance, state variable constraints 
(which may be based on physical considerations) are often neglected because they 
do not fit easily into the structure of the Kalman filter. This paper develops two an- 
alytic methods of incorporating state variable inequality constraints in the Kalman 
filter. The first method is a general technique of using hard constraints to enforce 
inequalities on the state variable estimates. The resultant filter is a combination of a 
standard Kalman filter and a quadratic programming problem. The second method 
uses soft constraints to estimate state variables that are known to vary slowly with 
time. (Soft constraints are constraints that are required to be approximately satis- 
fied rather than exactly satisfied.) The incorporation of state variable constraints 
increases the computational effort of the filter but significantly improves its estima- 
tion accuracy. The improvement is proven theoretically and shown via simulation 

* Corresponding author. This work was supported in part by a NASA/ASEE Summer 
Faculty Fellowship. 
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results. The use of the algorithm is demonstrated on a linearized simulation of a 
turbofan engine to estimate health parameters. The turbofan engine model con- 
tains 16 state variables, 12 measurements, and 8 component health parameters. It 
is shown that the new algorithms provide improved performance in this example 
over unconstrained Kalman filtering. 

1 Introduction 

For linear dynamic systems with white process and measurement noise, the Kalman 
filter is known to be an optimal estimator. However, in the application of Kalman 
filters there is often known model or signal information that is either ignored or 
dealt with heuristically [1], This paper presents two ways to generalize the Kalman 
filter in such a way that known inequality constraints among the state variables are 
satisfied by the state variable estimates. 

The first method presented here for enforcing inequality constraints on the state 
variable estimates uses hard constraints. It is based on a generalization of the ap- 
proach presented in [2] , which dealt, with the incorporation of state variable equality 
constraints in the Kalman filter. Inequality constraints are inherently more compli- 
cated than equality constraints, but standard quadratic programming results can be 
used to solve the Kalman filter problem with inequality constraints. At each time 
step of the constrained Kalman filter, we solve a quadratic programming problem 
to obtain the constrained state estimate. A family of constrained state estimates is 
obtained, where the weighting matrix of the quadratic programming problem deter- 
mines which family member forms the desired solution. It is stated in this paper, 
on the basis of [2], that the constrained estimate has several important properties. 
The constrained state estimate is unbiased and has a smaller error covariance than 
the unconstrained estimate. We show which member of all possible constrained so- 
lutions has the smallest error covariance. We also show the one particular member 
that is always (i.e., at each time step) closer to the true state than the unconstrained 
estimate. 
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The second method for enforcing inequality constraints uses soft constraints 
via a penalty term in the Kalman filter optimization problem. This prevents the 
state estimate from changing too rapidly. It essentially smooths the unconstrained 
Kalman filter estimate when the state variables are known to vary slowly with time. 
It is shown that the constrained state estimate is unbiased, approaches the uncon- 
strained estimate as time approaches infinity, and (under certain special conditions) 
is equal to the running average of the unconstrained estimate. 

The application considered in this paper is turbofan engine health parameter 
estimation [3]. The performance of gas turbine engines deteriorates over time. This 
deterioration reduces the fuel economy of the engine. Airlines periodically collect 
engine data in order to evaluate the health of the engine and its components. The 
health evaluation is then used to determine maintenance schedules. Reliable health 
evaluations are used to anticipate future maintenance needs. This offers the benefits 
of improved safety and reduced operating costs. The money-saving potential of such 
health evaluations is substantial, but only if the evaluations are reliable. The data 
used to perform health evaluations are typically collected during flight and later 
transferred to ground-based computers for post-flight analysis. Data are collected 
each flight at the same engine operating points and corrected to account for vari- 
ability in ambient conditions. Typically, data are collected for a period of about 
3 seconds at a rate of about 10 or 20 Hz. Various algorithms have been proposed 
to estimate engine health parameters, such as weighted least squares [4], expert 
systems [5], Kalman filters [6], neural networks [6], and genetic algorithms [7]. 

This paper applies constrained Kalman filtering to estimate engine component 
efficiencies and flow capacities, which are referred to as health parameters. We can 
use our knowledge of the physics of the turbofan engine in order to obtain a dynamic 
model [8, 9]. The health parameters that we try to estimate can be modelled as 
slowly varying biases. The state vector of the dynamic model is augmented to include 
the health parameters, which are then estimated with a Kalman filter [10]. The 
model formulation in this paper is similar to previous NASA work [11] . However, [11] 
was limited to a 3-state dynamic model and 2 health parameters, whereas this 
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present work includes a more complete 16-state model and 8 health parameters. In 
addition, we have some a priori knowledge of the engine’s health parameters: we 
know that they never improve. Engine health always degrades over time, and we can 
incorporate this information into state constraints to improve our health parameter 
estimation. (This is assuming that no maintenance or engine overhaul is performed.) 
This is similar to the probabilistic approach to turbofan prognostics proposed in [12]. 
It should be emphasized that in this paper we are confining the problem to the 
estimation of engine health parameters in the presence of degradation only. There 
are specific engine cases that can result in abrupt shifts in filter estimates, possibly 
even indicating an apparent improvement in some engine components. An actual 
engine performance monitoring system would need to include additional logic to 
detect and isolate such faults. 

Section 2 presents a discussion of the standard discrete time Kalman filter. Some 
important properties of the Kalman filter that will be used later in this paper are 
also reviewed. Section 3 generalizes the results of [2] to hard inequality constraints. 
This inequality-constrained Kalman filter has several attractive theoretical proper- 
ties, including state variable estimates that are unbiased, an estimation error vari- 
ance smaller than the unconstrained filter, and a time-domain estimation error that 
is always smaller than the unconstrained estimation error. Section 4 extends the 
standard Kalman filter in a different way for those cases where it is known that the 
state variables change slowly with time. This constraint is enforced by finding a new 
state estimate that is “close” to the unconstrained estimate in some sense, but that 
is slowly time varying. It is shown that this new estimate is unbiased, approaches 
the unconstrained estimate as time goes to infinity, and (under certain conditions) 
is equal to the running average of the unconstrained estimate. 

Section 5 discusses the problem of turbofan health parameter estimation, along 
with the dynamic model that we used in our simulation experiments. Although the 
health parameters are not state variables of the model, it is shown how the dynamic 
model can be augmented in such a way that a Kalman filter can estimate the health 
parameters [10, 11]. We then show how this problem can be expressed in such 
a way to be compatible with the constraints discussed in the preceding sections. 
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Section 6 presents some simulation results based on a turbofan model linearized 
around a known operating point. We show that the Kalman filter can estimate 
health parameters with an average error of less than 0.2%, and the constrained 
Kalman filters perform better than the unconstrained filter. Section 7 presents 
some concluding remarks and suggestions for further work. 

2 Kalman Filtering 

This section reviews standard (unconstrained) state estimation via the Kalman filter 
and some important properties of the filter that will be used later in this paper. The 
results and notation are taken from [13]. Consider the discrete linear time-invariant 
system given by 

Xk+i = Axk + Buk + Wk ( 1 ) 

Vk = Cx k + e k 

where k is the time index, x is the state vector, u is the known control input, y 
is the measurement, and {«>*.} and {e k } are noise input sequences. The problem 


is to find an estimate x k+ \ of x k+ i given the measurements {j/o> Vi, ■ ■ ■ , Vk}- We 
will use the symbol Y k to denote the column vector that contains the measurements 
{yo, yi, ■ ■ ■ ,y k }- We assume that the following standard conditions are satisifed. 

E[xo] = x 0 (2) 

E[w k ] = E[e k ] = 0 (3) 

E[(x 0 - x 0 )(x 0 - x 0 ) T ] = S 0 (4) 

E[WkWm] = Q S km (5) 

£[efcel] = RSkm (6) 

E[wu^n ] = E[x k e = E[x kW T\ = 0 (7) 


where E[-] is the expectation operator, x is the expected value of x, and 6 km is the 
Kronecker delta function ( 6 km = 1 if k = m, 0 otherwise). Q and R are positive 
semidefinite covariance matrices. The Kalman filter equations are given by 

K k = AE k C T (CE k C T + R)- 1 (8) 
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1 

= Ax k + Bu k + K k {y k - Cx k ) 

(9) 

Efc+i 

= {AY k -K k CY k )A T + Q 

(10) 


where the filter is initialized with xq — xq, and So given above. It can be shown [13] 
that the Kalman filter has several attractive properties. For instance, if Xq, {w k }, 
and {efc } are jointly gaussian, the Kalman filter estimate x k+ \ is the conditional 
mean of x k +\ given the measurements Yf; i.e., x k +i — E[x k +\\Y k \. Even if x$, 
{w k }, and {e k } are not jointly gaussian, the Kalman filter estimate is the best affine 
estimator given the measurements Yf; i.e., of all estimates of x k+ i that are of the 
form FY k + g (where F is a fixed matrix and g is a fixed vector) , the Kalman filter 
estimate is the one that minimizes the variance of the estimation error. It can be 
shown [13, pp. 92 ff.] that the Kalman filter estimate (i.e., the minimum variance 
estimate) can be given by 

®fc+l ®fc+l = -£fc+l “f '^-‘xyEyy (Y k Ff.) (11) 


where x k +\ is the mean of Xk+i, Y xy is the variance matrix of x k +\ and Y k , Y yy 
is the covariance matrix of Yf , and x k+ \ is the conditional mean of x k+ \ given 
the measurements Yf . In addition, from [13, p. 93] we know that the Kalman 
filter estimate x k+ \ and Yf are jointly gaussian, in which case x k+i is conditionally 
gaussian given Yf. The conditional probability density function of x k +i given Yf is 


_ exp[— (x — xj 1 £ l (.r-;f)/2] 
1 1 ’ (2-7r) n / 2 |E| 1 / 2 

where n is the dimension of x and 


( 12 ) 


E — Y xx — 'Bxy'Eyy Ti yx (13) 

The Kalman filter estimate is that value of x that maximizes the conditional prob- 
ability density function P{x\Y), and E is the covariance of the Kalman filter esti- 
mation error. 
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3 Kalman Filtering with Hard Inequality Con- 
straints 

This section extends the well known results of the previous section to cases where 
there are known linear inequality constraints among the state components. Also, 
several important properties of the constrained filter are discussed. Consider the 
dynamic system of (1) where we are given the additional constraint 

Dx k < d k (14) 

where D is a known s x n constant matrix, s is the number of constraints, n is the 
number of state variables, and s < n. It is assumed in this paper that D is full 
rank, i.e. , that D has rank s. This is an easily satisfied assumption. If D is not full 
rank that means we have redundant state constraints. In that case we can simply 
remove linearly dependent rows from D (i.e., remove redundant state constraints) 
until D is full rank. Three different, approaches to the constrained state estimation 
problem are given in this section. The time index k is omitted in the remainder of 
this section for ease of notation. 

3.1 The Maximum Probability Method 

In this section we derive the constrained Kalman filtering problem by using a max- 
imum probability method. From [13, pp. 93 ff.] we know that the Kalman filter 
estimate is that value of x that maximizes the conditional probability density func- 
tion P(x\Y), which is given in (12). The constrained Kalman filter can be derived 
by finding an estimate x such that the conditional probability P(x\Y) is maximized 
and x satisfies the constraint (14). Maximizing P{x\Y) is the same as maximizing 
its natural logarithm. So the problem we want to solve can be given by 

maxlnP(5|K) =$■ min(5 — x) t YT 1 (x — x) (15) 

X 

such that Dx < d 
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Using the fact that the unconstrained state estimate x — x (the conditional mean 
of x), we rewrite the above equation as 

min(5 T S _1 5 — 2x T ’E~ 1 x) such that Dx < d (16) 

X 

Note that this problem statement depends on the conditional gaussian nature of x, 
which in turn depends on the gaussian nature of xq, {it;*.}, and {e*,} in (1). 

3.2 The Mean Square Method 

In this section we derive the constrained Kalman filtering problem by using a mean 
square minimization method. We seek to minimize the conditional mean square 
error subject to the state constraints. 

min E(\\x — x |j 2 |K) such that Dx < d (17) 

X 

where || • | denotes the vector two-norm. If we assume that x and Y are jointly 
gaussian, the mean square error can be written as 

E{\\x — x || 2 |F) = J (x — x) 1 (x — x)P(x\Y)dx (18) 

= J x T xP(x\Y)dx — 2x J xP{x\Y)dx + x T x (19) 

Noting that the Kalman filter estimate is the conditional mean of x, i.e., 

x = j xP{x\Y)dx (20) 

we formulate the first order conditions necessary for a minimum as 

vYm(x T x — 2x T x) such that Dx < d (21) 

X 

Again, this problem statement depends on the conditional gaussian nature of x, 
which in turn depends on the gaussian nature of xq, {iCfc}, and {e*,} in (1). 

3.3 The Projection Method 

In this section we derive the constrained Kalman filtering problem by directly pro- 
jecting the unconstrained state estimate x onto the constraint surface. That is, we 
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solve the problem 


min(5 — x) T W ( x — x ) such that Dx < d (22) 

X 

where W is any symmetric positive definite weighting matrix. This problem can be 
rewritten as 

Tam{x T Wx — 2x7 IV x) such that Dx < d (23) 

X 

The constrained estimation problems derived by the maximum probability method (16) 
and the mean square method (21) can be obtained from this equation by setting 
W — E -1 and W — I respectively. Note that this derivation of the constrained 
estimation problem does not depend on the conditional gaussian nature of x; i.e., 
x 0 , {wk}, and {e*,} in (1) are not assumed to be gaussian. 

3.4 The Solution of the Constrained State Estimation 
Problem 

The problem defined by (23) is known as a quadratic programming problem [14, 15]. 
There are many algorithms for solving quadratic programming problems, almost all 
of which fall in the category known as active set methods. An active set method 
uses the fact that it is only those constraints that are active at the solution of the 
problem that are significant in the optimality conditions. Assume that t of the s 
inequality constraints are active at the solution of (23) , and denote by D and d the t 
rows of D and t elements of d corresponding to the active constraints. If the correct 
set of active constraints was known a priori then the solution of (23) would also be 
a solution of the equality-constrained problem 

min(5 T IT. ; r — 2x T Wx) such that Dx = d (24) 

X 

This shows that the inequality constrained problem defined by (23) is equivalent to 
the equality-constrained problem defined by (24). The equality-constrained problem 
was discussed in [2], and so those results can be used to investigate the properties 
of the inequality-constrained problem. 
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3.5 Properties of the Constrained State Estimate 

In this section we examine some of the statistical properties of the constrained 
Kalman filter. We use x to denote the state estimate of the unconstrained Kalman 
filter, and x to denote the state estimate of the constrained Kalman filter as given 
by (23), recalling that (16) and (21) are special cases of (23). 

Theorem 1 The solution x of the constrained state estimation •problem given by (23) 
is an unbiased state estimator for the system (1) for any symmetric positive definite 
weighting matrix W . That is, 

E{x) — E{x ) (25) 

Theorem 2 The solution x of the constrained state estimation problem given by (23) 
with W = S _1 , where S is the covariance of the unconstrained estimate given in (10) 
and (13), has an error covariance that is less than or equal to that of the uncon- 
strained state estimate. That is, 

Cov{x — x) < Cov(x — x) (26) 


At first this seems counterintuitive, since the standard Kalman filter is by definition 
the minimum variance filter. However, we have changed the problem by introducing 
state variable constraints. Therefore, the standard Kalman filter is no longer the 
minimum variance filter, and we can do better with the constrained Kalman filter. 

Theorem 3 Among all the constrained Kalman filters resulting from the solution 
of (23), the filter that uses W = E _1 has the smallest estimation error covariance. 
That is, 

Cov(x%- 1 ) < Cov(x w ) for all W (27) 

Theorem 4 The solution x of the constrained state estimation problem given by (23) 
with W = I satisfies the inequality 

II X k - x k \\ < \\x k - x k \\ for all k (28) 

where || ■ | is the vector two-norm and x is the unconstrained Kalman filter estimate. 
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Theorem 5 The error of the solution x of the constrained state estimation problem 
given by (23) with W — I is smaller than the unconstrained estimation error in the 
sense that 

Tr[Cov(x)\ < Tr[Cov(x)\ (29) 

where Tr[-] indicates the trace of a matrix, and Cov(-) indicates the covariance matrix 
of a random vector. 

The above theorems all follow from the equivalence of (23) and (24), and the 
proofs presented in [2], We note that if any of the s constraints are active at the 
solution of (23), then strict inequalities hold in the statements of Theorems 2-5. The 
only time that equalities hold in the theorems is if there are no active constraints at 
the solution of (23); that is, if the unconstrained Kalman filter satisfies the inequality 
constraints. 

4 Kalman Filtering with Soft Inequality Con- 
straints 

In this section we are interested in obtaining a Kalman filter-based state estimate 
for state variables which we know a priori vary slowly with time. Since we are 
concerned with using the Kalman filter as a parameter estimator, we will assume for 
this problem that the A matrix in (1) is the identity matrix and the B matrix is zero. 
With this in mind, we can use the results of the previous section, especially (22), to 
formulate a Kalman filter-based estimate as follows 

min(.Tfc — x k ) T W(x k ~ %k) such that x{i } varies slowly (30) 

•Efc 

where, as before, W is a constant symmetric positive definite weighting matrix. This 
is a type of regularization; that is, some additional structure is incorporated into 
the Kalman filter estimate [16, 17, 18]. The above problem can be formulated as 

min[(.T fc - x k ) T W(x k - x k ) + ( x k - x k _i) T V k (x k - £ fc -i)] (31) 

k 
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where V*. is a (possibly time-varying) symmetric positive definite weighting matrix 
that balances the desire for a close approximation to x and smooth estimate x. The 
solution to the above problem is 

x Q = E[x o] (32) 

x k = (I V + V k )-\Wx k +V k x k _ 1 ) 

Since W and V k are both positive definite, we know that (W + k/,-) -1 exists. 

Theorem 6 Assume (as stated above) that A — I and B — 0 in (1). Then the 
solution x of the constrained state estimation problem given by (32) is an unbiased 
state estimator for the system (1) for any symmetric positive definite weighting 
matrices W and V k . That is, 

E(x) — E(x ) (33) 

Proof: The theorem can be proven by induction. Since A — I and B — 0 we know 
that E[x k \ — xq for all k. We therefore know from (32) that xq = xq. From (32) 
with k = 0 we see that E[x\] — xq. We repeat this process to show that E[x k \ — 
E[x k \ = xq for all k. 

QED 

Theorem 7 Assume (as stated above) that A — I and B — 0 in (1). Further 
assume that w k = 0 in (1) (since we are trying to estimate constant parameters). 
Then the constrained state estimate x approaches the unconstrained estimate x in 
the limit as time goes to infinity. That is, 

lim x k = lim x k (34) 

k — »oo k — >oo 

Proof: We see from (8)-(10) that, under the conditions stated here, K k — > 0 as 
k — > oo. Therefore x k approaches a constant value as k — > oo. From (32) we see 
that, in steady state 

x = (W + V k )~ 1 (Wx + V k x) (35) 

= [I-iW+V^V^iW + V^Wx 
= {I + W~ l V k )(W + VkT'Wx 
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where the last equality follows from the matrix inversion lemma. Premultiplying 
both sides of the above equation by W we obtain Wx — Wx, so if W is invertible 
(which it is, since we are assuming in this section that W is positive definite), we 
obtain x — x (in steady state). Note that the theorem is true even if Vk does not 
approach a steady state value as k — > oo. 

QED 

Theorem 8 If Vk — {k — 1 )W in (32) then Xk is the running average of Xk- 

Proof: The running average of Xk is defined as 

k 

Xk = \±* (36) 

i= 1 

which implies that 

Xk+i — ^ j (®fc+i + kXk ) (37) 

Now if Vk = (k — l)W then (32) shows that 

x k+ i = [(k + l)W]-\Wxk+i+kWx k ) (38) 

= Y^{x k +i + kx k ) 

which is exactly the running average shown in (37). 

QED 

5 Turbofan Engine Health Monitoring 

Figure 1 shows a schematic representation of a turbofan engine. A single inlet 
supplies airflow to the fan. Air leaving the fan separates into two streams: one 
stream passes through the engine core, and the other stream passes through the 
annular bypass duct. The fan is driven by the low pressure turbine. The air passing 
through the engine core moves through the compressor, which is driven by the high 
pressure turbine. Fuel is injected in the main combustor and burned to produce 
hot gas for driving the turbines. The two air streams combine in the augmentor 
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duct, where additional fuel is added to further increase the air temperature. The air 
leaves the augmentor through the nozzle, which has a variable cross section area. 

Various turbofan simulation packages have been proposed over the years [19, 20, 
21]. This model is based on a gas turbine engine simulation software package called 
DIGTEM (Digital Turbofan Engine Model) [8, 22]. DIGTEM is written in Fortran 
and includes 16 state variables. It uses a backward difference integration scheme 
because the turbofan model contains time constants that differ by up to four orders 
of magnitude. 

The nonlinear equations used in DIGTEM can be found in [8, 9]. The time- 
invariant equations can be summarized as follows. 

x = f(x,u,p) + wi(t) (39) 

V = g(x,u,p)+e(t) 

x is the 16-element state vector, u is the 6-element control vector, p is the 8-element 
vector of health parameters, and y is the 12-element vector of measurements. The 
noise term w\ (t) represents inaccuracies in the model, and e(t) represents measure- 
ment noise. The elements in these vectors are summarized in Tables 1-4, along with 
their values at the nominal operating point (xo,tto,Po,yo) considered in this paper. 
Table 4 also shows typical signal- to- noise ratios for the measurements, based on 
NASA experience and previously published data [23]. Sensor dynamics are assumed 
to be high enough bandwidth that they can be ignored in the dynamic equations [23] . 
Equation (39) can be linearized about the nominal operating point by using the first 
order approximation of the Taylor series expansion 

f(x,u,p) « f(x 0 ,u 0 ,p 0 )+ (40) 

g(x,u,p) ■ ' g(x 0 ,uo,po) + ~^r( x ~ x o) + ~g~( u ~ °o) + -]4kp-Po) + 

Therefore, a linear small signal system model can be defined for small excursions 
from the nominal operating point. 
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8x = x — x o 

fiy = y-yo 


AiSx + B8u + A-26p + wi(t) 
C\8x + D8u + C-28p + e{t) 


(41) 


We note that 


A, 

M (i,j) 


d£ 

dx 

A x(i) 
Ax(j) 


(42) 


Similar equations hold for the A-2, C\, and C2 matrices. We obtained numerical 
approximations to the A\, A2, C\, and C 2 matrices by varying x and p from their 
nominal values (one element at a time) and recording the new x and y vectors in 
DIGTEM. 

Turbofan engine health monitoring is typically a two-step process [3]. In the 
first step, engine data is collected each flight at the same engine operating points 
and corrected to account for variability in ambient conditions. Data are typically 
collected for a period of about 3 seconds per flight at a rate of about 10 or 20 Hz. In 
the second step, the data are transferred to ground-based computers for post-flight 
analysis to determine engine health. 

The goal of our turbofan engine health monitoring problem is to obtain an 
accurate estimate of 8p, which varies slowly with time. We therefore assume that 
8p is constant between measurement times. We also assume that the control input 
is perfectly known, so 8u = 0. This gives us the following equivalent discrete time 
system [24, pp. 90 ff.]. 


8x k+1 = A ld 8x k + A 2d 8pk + w lk (43) 

8yk = C\8x k + C 2 8p k + e k 

where Ai d = exp(AiT) and A 2d = A ^ 1 (A\ d — I)A2 (assuming that A\ is invertible, 
which it is in our problem). We next augment the state vector with the health 
parameter vector [11] to obtain the system equation 
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(44) 


fo'Ek+l 

fiPk+1 

Sy k 


Aid 

A-2d 


6x k 

+ 

W\k 

0 

I 


Sp k 

W'2k 

[<* 

c 2 } 

6x k , „ 
6* +e ‘ 



where w 2 k is a small noise term (uncorrelated with W\ k ) that represents model 
uncertainty and allows the Kalman filter to estimate time-varying health parameter 
variations. The discrete time small signal model can be written as 


7 r k + 1 
fiPk+l 

% 



c 


Sx k 

Spk 


+ Wk 
+ e k 


(45) 


where the definitions of A and C are apparent from a comparison of the two preced- 
ing equations. Now we can use a Kalman filter to estimate 6x k and 6p k . Actually, 
we are only interested in estimating 6p k (the health parameter deviations), but the 
Kalman filter gives us the bonus of also estimating 6x k (the excursions of the original 
turbofan state variables). 

It is known that health parameters do not improve over time. That is, 6p( 1), 
6p( 2), 6p( 3), 6p{ 4), 6p{ 6), and 6p( 8) are always less than or equal to zero and 
always decrease with time. Similarly, 6p{ 5) and Sp(7) are always greater than or 
equal to zero and always increase with time. In addition, it is known that the health 
parameters vary slowly with time. As an example, since 6p( 1) is the constrained 
estimate of <5p(l), we can enforce the following constraints on 6p( 1). 


6p( 1) < 0 (46) 

fy>k+ 1(!) < %(!) + 7i + 

fy>k+ 1(!) > %(l)-7f 


where qj 1 " and qf are nonnegative factors chosen by the user that allows the state 
estimate to vary only within prescribed limits. Typically we choose qj“ > qj* - so that 
the state estimate can change more in the negative direction than in the positive 
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direction. This is in keeping with our a priori knowledge that this particular state 
variable never increases with time. Ideally we would have 7 + = 0 since <5p(l) never 
increases. However, since the state variable estimate varies around the true value of 
the state variable, we choose 7 ^" > 0. This allows some time-varying increase in the 
state variable estimate to compensate for a state variable estimate that is smaller 
than the true state variable value. 

These constraints are linear and can therefore easily be incorporated into the 
form required in the constrained filtering problem statement (14). Note that this 
does not take into account the possibility of abrupt changes in health parameters 
due to discrete damage events. That possibility must be addressed by some other 
means (e.g., residual checking [3]) in conjuction with the methods presented in this 
paper. 


6 Simulation Results 

We simulated the methods discussed in this paper using MATLAB. We simulated 
a steady state 3 second burst of engine data measured at 10 Hz during each flight. 
Each of these routine services was performed at the single operating point shown 
in Tables 1-4. The signal-to-noise ratios were determined on the basis of NASA 
experience and previously published data [23] and are shown in Table 4. We used a 
one-sigma process noise in the Kalman filter equal to 1% of the nominal state values 
to allow the filter to be responsive to changes in the state variables. We set the one 
sigma process noise for each component of the health parameter portion of the state 
derivative equation to 0.01% of the nominal parameter value. This was obtained by 
tuning. It was small enough to give reasonably smooth estimates, and large enough 
to allow the filter to track slowly time-varying parameters. For the filter with hard 
constraints, we chose the 7 variables in (46) such that the maximum allowable rate 
of change in 6p was a linear 9% per 500 flights in the direction of expected change, 
and 3% per 500 flights in the opposite direction. The true health parameter values 
never change in a direction opposite to the expected change. However, we allow 
the state estimate to change in the opposite direction to allow the Kalman filter 
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to compensate for the fact that the state estimate might be either too large or too 
small. We set the weighting matrix W in (23) and (31) equal to E _1 in accordance 
with Theorem 3. We found by experimenting that setting the weighting matrix Vk 
in (31) equal to 120 W resulted in good performance for the Kalman filter with soft 
constraints. 

The first test case we simulated was a linear degradation of the first health 
parameter (fan airflow) over 500 flights, while the other seven health parameters 
remained constant. Figure 2 shows the Kalman filters’ performances in this case. We 
ran eight simulations like this. In each simulation, one of the eight health parameters 
degraded linearly by a factor of 3% during the course of the simulation, while the 
other seven health parameters remained constant. The 3% degradation over 500 
flights is in line with turbofan performance data collected by NASA and reported in 
the literature [25]. Each of the eight cases exhibit performance similar to Figure 2. 
Table 5 shows the performance of the filters averaged over all eight simulations. 
All of the filters estimate the health parameters to within less than 0.2% of their 
nominal values. It can be seen that (on average) the filter with soft constraints 
offers an 11% improvement over the unconstrained filter, and the filter with hard 
constraints offers a 22% improvement over the unconstrained filter. These numbers 
should not be interpreted as having any statistical signficance (due to our limited 
sample size of eight cases) but they do show the improvement that is possible with 
constrained Kalman filters. Table 5 also shows that a couple of health parameters 
(fan airflow and LPT airflow) were actually estimated better with the unconstrained 
filter than with the constrained filter. We therefore see that the constrained filter 
does not guarantee better estimation in every individual sample run, but it does 
guarantee better performance statistically. 

The next scenario we considered was the case where all eight health param- 
eters degrade at the same time. We simulated a degradation over 500 flights of 
— 1% for fan airflow, —2% for fan efficiency, —3% for compressor airflow, —2% for 
compressor efficiency, +3% for high pressure turbine airflow, —2% for high pressure 
turbine enthalpy change, +2% for low pressure turbine airflow, and —1% for low 
pressure turbine enthalpy change. This is summarized in Table 6. Figure 3 shows 
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the performance of the Kalman filters in this case. Table 7 shows the performance 
of the filters averaged over 16 simulations like this (each simulation being subject 
to a different random noise history). It can be seen that (on average) the filter with 
soft constraints offers a 9% improvement over the unconstrained filter, and the filter 
with hard constraints offers a 38% improvement over the unconstrained filter. As 
mentioned above, these numbers should not be interpreted as having any statisti- 
cal signficance (due to our limited sample size of 16 cases) but they do show the 
improvement that is possible with constrained Kalman filters. 

The improved performance of the constrained filters comes with a price, and 
that price is computational effort. The filter with soft constraints requires only 
slightly (14%) more computational effort than the unconstrained filter, but the 
filter with hard constaints requires about four times the computational effort of 
the unconstrained filter. This is because of the additional quadratic programming 
problem that is required for hard constraints. However, computational effort is not 
a critical issue for the particular application of turbofan health estimation since the 
filtering is performed on ground-based computers after each flight. 

7 Conclusion and Discussion 

We have presented two methods for incorporating linear state inequality constraints 
in a Kalman filter. The first method incorporated hard constraints into the Kalman 
filter to maintain the state variable estimates within a user-defined envelope. The 
second method incorporated soft constraints into the Kalman filter to ensure that 
the state variable estimates vary slowly with time. The simulation results demon- 
strate the effectiveness of these methods, particularly for turbofan engine health 
estimation. 

If the system whose state variables are being estimated has known state variable 
constraints, then those constraints can be incorporated into the Kalman filter as 
shown in this paper. However, in practice, the constraints enforced in the filter 
might be more relaxed than the true constraints. This allows the filter to correct 
state variable estimates in a direction that the true state variables might never 
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change. This is a departure from strict adherence to theory, but in practice this 
improves the performance of the filter. This is an implementation issue that is 
conceptually similar to tuning a standard Kalman filter. 

It was seen in Theorem 2 that the filter with hard constraints has a smaller 
estimation error covariance than the unconstrained Kalman filter. At first this seems 
counterintuitive, since the standard Kalman filter is by definition the minimum 
variance filter. However, we have changed the problem by introducing state variable 
constraints. Therefore, the standard Kalman filter is not the minimum variance 
filter for the turbofan engine health estimation problem, and we can do better with 
the constrained Kalman filter. 

We saw that the filter with hard constraints required a much larger computa- 
tional effort than the standard Kalman filter. This is due to the addition of the 
quadratic programming problem that must be solved in the constrained Kalman 
filter. The engineer must therefore perform a tradeoff between computational ef- 
fort and estimation accuracy. For real time applications the improved estimation 
accuracy may not be worth the increase in computational effort. 

It was seen in Figures 2 and 3 that although the constrained filters improve 
the estimation accuracy, the general trend of the state variable estimates does not 
change with the introduction of state constraints. This is because the constrained 
filters are based on the unconstrained Kalman filter. The constrained filter estimates 
therefore have the same shape as the unconstrained estimates until the constraints 
are violated, at which point the state variable estimates are projected onto the 
edge of the constraint boundary. The constrained filters presented in this paper 
are not qualitatively different than the standard Kalman filter; they are rather a 
quantitative improvement in the standard Kalman filter. 

Note that the Kalman filter works well only if the assumed system model matches 
reality fairly closely. The method presented in this paper, by itself, will not work 
well if there are large sensor biases or hard faults due to severe component failures. 
A mission-critical implementation of a Kalman filter should always include some sort 
of residual check to verify the validity of the Kalman filter results, particularly for 
the application of turbofan engine health estimation considered in this paper [3, 26]. 


NASA/TM— 2003-212111 


20 



Although we have considered only linear state constraints, it is not conceptually 
difficult to extend this paper to nonlinear constraints. If the state constraints are 
nonlinear they can be linearized as discussed in [2]. 

Further work along the lines of this research could focus on combining our work 
with [27] in order to guarantee convergence in the presence of nonlinear constraints. 
Other efforts could explore the incorporation of state constraints for optimal smooth- 
ing, or the use of state constraints in H 0 0 filtering [28]. Further work could also focus 
on integrating the nonlinear simulation logic in DIGTEM [8, 22] with the Kalman 
filter to obtain more complete results. This would also allow us to more easily 
test the Kalman filter at various operating points without translating data from 
DIGTEM to MATLAB. 
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State 

Nominal Value 

Low Pressure Turbine Rotor Speed 

6140 RPM 

High Pressure Turbine Rotor Speed 

9395 RPM 

Compressor Mass Flow 

0.457 kg/s 

Combustor Inlet Temperature 

965 K 

Combustor Mass Flow 

0.264 kg/s 

High Pressure Turbine Inlet Temperature 

1593 K 

High Pressure Turbine Mass Flow 

1.48 kg/s 

Low Pressure Turbine Inlet Temperature 

1129 K 

Low Pressure Turbine Mass Flow 

1.79 kg/s 

Augmentor Inlet Temperature 

790 K 

Augmentor Mass Flow 

1.46 kg/s 

Nozzle Inlet Temperature 

790 K 

Duct Fluid Momentum 

53.6 kg/s 2 

Augmentor Fluid Momentum 

103 kg/s 2 

Duct Mass Flow 

4.52 kg/s 

Duct Temperature 

571 K 


Table 1: Turbofan states. 


Control 

Nominal Value 

Combustor Fuel Flow 
Augmentor Fuel Flow 
Nozzle Throat Area 
Nozzle Exit Area 
Fan Vane Angle 
Compressor Van Angle 

0.37 kg/s 
0 kg/s 
430 cm 2 
492 cm 2 
-25 deg 
-20 deg 


Table 2: Turbofan controls. 
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Health Parameter 

Nominal Value 

Fan Airflow 

Fan Efficiency 

Compressor Airflow 

Compressor Efficiency 

High Pressure Turbine Airflow 

High Pressure Turbine Enthalpy Change 

Low Pressure Turbine Airflow 

Low Pressure Turbine Enthalpy Change 

102 kg/s 
0.82 

48.7 kg/s 
0.83 

41.0 kg/s 
101 J/kg 
48.3 kg/s 

27.1 J/kg 


Table 3: Turbofan health parameters. 


Measurement 

Nominal Value 

SNR 

Low Pressure Turbine Rotor Speed 

6140 RPM 


High Pressure Turbine Rotor Speed 

9395 RPM 


Duct Pressure 

19.0 N/cm 2 

200 

Duct Temperature 

571 K 

100 

Compressor Inlet Pressure 

20.5 N/cm 2 


Compressor Inlet Temperature 

577 K 

100 

Combustor Pressure 

97.5 N/cm 2 

200 

Combustor Inlet Temperature 

965 K 

100 

Low Pressure Turbine Inlet Pressure 

26.8 N/cm 2 

100 

Low Pressure Turbine Inlet Temperature 

1130 K 

70 

Augment or Inlet Pressure 

17.4 N/cm 2 

100 

Augmentor Inlet Temperature 

790 K 

70 


Table 4: Turbofan measurements. 
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Estimation Error (%) 

Health 

Unconstrained 

Soft Constrained 

Hard Constrained 

Parameter 

Filter 

Filter 

Filter 

Fan Airflow 

0.123 


0.139 

Fan Efficiency 

0.177 


0.113 

Compressor Airflow 

0.145 


0.113 

Compressor Efficiency 

0.102 

0.086 

0.059 

HPT Airflow 

0.116 

0.100 

0.101 

HPT Enthalpy Change 

0.093 

0.081 

0.055 

LPT Airflow 

0.104 

0.090 

0.109 

LPT Enthalpy Change 

0.181 

0.168 

0.118 

Average 

0.130 

0.116 

0.101 


Table 5: Kalman filter estimation errors. HPT = High Pressure Turbine, 
and LPT = Low Pressure Turbine. The numbers shown are RMS estimation 
errors (percent) averaged over eight simulations where each simulation had one 
health parameter degradation while the other seven health parameters were 
unchanged. 


Health 

True Degradation 

Fan Airflow 

-1% 

Fan Efficiency 

-2% 

Compressor Airflow 

-3% 

Compressor Efficiency 

- 2 % 

HPT Airflow 

+3% 

HPT Enthalpy Change 

-2% 

LPT Airflow 

+2% 

LPT Enthalpy Change 

-1% 


Table 6: Health parameter degradation amounts for test scenario. 
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Estimation Error (%) 


Health 

Parameter 

Unconstrained 

Filter 

Soft Constrained 
Filter 

Hard Constrained 
Filter 

Fan Airflow 

0.129 

0.113 


Fan Efficiency 

0.163 

0.149 

0.105 

Compressor Airflow 

0.152 

0.146 

0.103 

Compressor Efficiency 

0.101 

0.087 

0.052 

LIPT Airflow 

0.119 

0.114 

0.076 

HPT Enthalpy Change 

0.092 

0.078 

0.050 

LPT Airflow 

0.104 

0.091 

0.057 

LPT Enthalpy Change 

0.168 

0.155 

0.111 

Average 

0.128 

0.116 

0.080 


Table 7: Kalman filter estimation errors. HPT = High Pressure Turbine, 
and LPT = Low Pressure Turbine. The numbers shown are RMS estimation 
errors (percent) averaged over 16 simulations, where each simulation had a 
linear degradation of all eight health parameters. 
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Figure 1 : Schematic representation of turbofan engine 
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(a) Unconstrained 
Kalman filter 



(b) Kalman filter 
with soft constraints 



500 


100 200 300 

flight number 


400 


e-0.5 


-2.5 


(c) Kalman filter with 
hard constraints 



Figure 2: Kalman filter estimates of health parameters. The fine health parameter changes 
were a -3% change in the first parameter, and zero change in the other seven parameters. 
The true health parameter changes are shown as heavy lines, and the filter estimates are 
shown as lighter lines. 
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(a) Unconstrained 
Kalman filter 



flight number 


(b) Kalman filter 
with soft constraints 



(c) Kalman filter with 
hard constraints 



Figure 3: Kalman filter estimates of health parameters. The fine health parameter changes 
were various values in between -3% and +3%. The fine health parameter changes are 
shown as heavy lines, and the filter estimates are shown as lighter lines. 
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